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An approach is presented for the atomistic study of phonon transport in real dielectric nanowires 
via Green functions. The formalism is applied to investigate the phonon flow through nanowires 
coated by an amorphous material. Examples for a simple model system, and for real Si nanowires 
coated by silica are given. New physical results emerge for these systems, regarding the character 
of the transition from ballistic to diffusive transport, the low temperature thermal conductance, 
and the influence of the wire-coating interface on the thermal transport. An efficient treatment of 
phonon scattering by the amorphous coating is also developed, representing a valuable tool for the 
investigation of thermal conduction through amorphous-coated nanowires. 
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I. INTRODUCTION 

The problem of phonon transport through dielectric 
wires of nanometer thickness^ is of special interest at 
present, since accurate measurements of thermal conduc- 
tion in these systems begin to be availableSii4. Theo- 
retically, this problem has been considered by different 
approaches. Some of these are the Boltzmann transport 
equation (BTE)^ Molecular Dynamics (MD)fi^i, and the 
transmission function approacl&SiiSiiiii^ii^. The BTE 
has been succesfuUy used in nanowire transport at high 
temperatures, for systems where the resistive length of 
the wire is long enough for the transport to be diffusive. 
Molecular Dynamics has the advantage that it can accu- 
rately consider the anharmonic interactions between the 
atoms. MD was applied to nanowires in ref. for exam- 
ple. It is nonetheless difficult to study low temperatures 
using MD, since it provides a classical description of the 
system. 

The transmission function approach is very well suited 
to study cases when phonons flow ballistically or semi- 
ballistically. In the limit of very low temperatures, an 
elastic continuum model provides a good description of 
the phonon flow. This model has yielded new insight into 
the problem of phonon scattering by surface roughness, 
for example, where the transition from ballistic to dif- 
fusive transport and the low temperature "dip" in the 
quantized thermal conductanceii have been investigated. 
Transmission function approaches have also been succes- 
fuUy applied to study the heat flow through a mesoscopic 
link^^, a nanocristal^"*, and monatomic chains^. Despite 
these excellent works, and unlike the case of electron 
transport, the use of the transmission function approach 
in phonon transport is still scarce. 

In this paper we develop a novel transmission function 
approach, and apply it to study a problem that has not 
been theoretically investigated before: the phonon trans- 
port along nanowires in which part of the length is sur- 
rounded by a thick coating of amorphous material. This 
is an important problem, since most dielectric nanowires 
are naturally coated by a layer of amorphous materiaU*i^. 



After explaining the general formalism, its use is il- 
lustrated by the study of phonon transmission through 
two concrete systems: a 1-D chain partially coated with 
an amorphous layer, and real 3-D nanowires with part 
of their length coated by a silica layer. Three main phe- 
nomena are investigated: 

-the transition from ballistic to diffussive transport. 
Explicit curves are shown of how this transition takes 
place, as a function of frequency and coated length. 

-the low temperature conductance. It is found that a 
phenomenon similar to the " dip" occuring in rough-edge 
uncoated wires^^, occurs also in amorphous coated wires, 
although displaying some differences with respect to the 
former case. 

-the effect of the wire-coating interface. We find an 
interesting saturation effect as a function of the coupling 
strength between the atoms at the interface. We also 
show the important influence of the interatomic link num- 
ber and configuration on the phonon transmission. 

Our approach is also new in many respects: 

1) An atomistic description is used. This is necessary 
since we are interested in how the atomic structure of 
the interface between the wire and the coating affects 
transport. This also enables to consider the whole dis- 
persive spectrum, and not only the lower frequencies. 
Although excellent atomistic investigation of transport 
through monatomic chains was done in ref. 0, we have 
not been able to find any other transmission function 
atomistic study of larger systems, like the Si nanowires 
considered here. Despite the method is similar to that of 
tight binding electron transport, it is important to pro- 
vide a complete and independent derivation of the formu- 
las in the phonon framework, specially for those readers 
not familiar with the theory of electron transport. Part 
of the paper is devoted to do this in detail. 

2) We present a new method that allows to treat the 
amorphous overlayer in an efficient manner. As we will 
see, the presence of the amorphous overlayer introduces 
a problem of overwhelming computational demand. We 
develop a technique that is able to yield the thermal con- 
ductance in the limit of thick overlayer, with much less 
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computational effort. The validity of this approach is 
explicitly shown with one example. 

The structure of the paper is as follows. Section ^ex- 
plains the formalism. An introduction to the transmis- 
sion function concept fScct. IlL^ is followed by the core 
of the method, Section FlI Bl where we derive the formu- 
las to compute the phonon transmission function. After- 
wards, Section III CI describes the specific self-consistent 
technique that enables us to treat the amorphous coat- 
ing surrounding the wires. Section lllll shows results for 
the two examples treated. Conclusions are summarized 
in Section Hvl 



a case we no longer have translational periodicity, and 
we cannot define anything like a "phonon band disper- 
sion" . Nevertheless, the concept of transmission chan- 
nels is still valid. As it was shown by Landauer— , no 
matter how scatterers are distributed in the system, one 
can always calculate a transmission fmiction, S(cj), which 
describes the propagation of quasiparticles between two 
reservoirs connected to the system, with different chemi- 
cal potentials. In contrast to translationally periodic sys- 
tems, where the conduction channels always yield integer 
quanta of transmission, now the waves can be partially 
reflected resulting in non integral transmissions. Thus in 
this case, instead of Eq. we have 



II. METHOD 

A. The transmission function in thermal 
conductivity and conductance 

Both the thermal conductivity and the thermal con- 
ductance can be obtained by the same single approach. 
Let us imagine a perfectly harmonic and translationally 
invariant nanowire, totally free of defects. The wire is 
free standing, so that there is no heat leakage from the 
boundary. The harmonicity assumption implies that, un- 
less there is scattering due to disordered defects or impu- 
rities, the thermal conductance of the whole wire would 
be independent of its length. In this situation, the heat 
flux associated with a small temperature difference AT 
between the two wire ends is given by the sum of the 
contributions of the individual phononic subbandsiS: 
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where a labels the phonon subbands, a is the wire's unit 
cell length, ks is Boltzmann's constant, k is the wave vec- 
tor in the axial direction, uJa{k) is the dispersion relation 

of subband a, u;™'*^''^*"^ is the lower (upper) frequency 
limit of the subband, and oj is the frequency. We note 
that the phonons inside the wire have a non equilibrium 
distribution different from the the Bose-Einstein type. 
As it is clear from the second line, it is not necessary to 
know the dependence of uja on k, but only the frequency 
limits of the subband. 

The (length independent) thermal conductance for this 
limit is then 
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In general, however, one has a nanowire with disor- 
dered scattering sources along a finite length. In such 
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where the phonon spectrum extends between frequencies 
and w/™. Eq. @ is of general validity and involves no 
approximations. (This equation is rigorously derived in 
Section III B 31 ) In the next subsection we explain how S 
is calculated exactly for independent phonons. 

Scattering causes the conductance to vary with the 
nanowire's length. For 3-D systems in the diffusive 
regime, the conductance varies inversely proportional to 
the system's length, L. Because of this one usually de- 
fines a " conductivity" , as 



Lcr(L) 



(3) 



where s is the sample's cross section. However, for 
nanowires, the length dependence of a cannot be just 
assumed to have a 1/L form, but needs to be calculated. 
For short L, a saturates to a sat- For long L the behavior 
depends on the type of scattering and the wire properties. 
The transition between the two limits also constitutes an 
important problem (studied in Section llll|) . In order to 
attack these issues, the transmission function has to be 
calculated for the atomically described, non-periodical, 
infinite system. The formalism allowing to do this is ex- 
plained in the next subsection. 



B. Calculation of the transmission function 

1. Interatomic potentials and the dynamical matrix 

The system's motion is determined by its dynamical 
matrix, obtained from the interatomic potentials of the 
system^. The nondiagonal elements of the dynamical 
matrix K are calculated as 



duiduj 



(4) 



where E is the energy, and Ui is the displacement of 
the i"* degree of freedom with respect to its equilibrium 
value. The diagonal elements are ka = X^j^^i —kij^^. 
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The dynamical equation of the system is 

{lu^M + K)u = (5) 

where M is a diagonal matrix with elements correspond- 
ing to the masses of the constituent atoms, and w is 
the vibrational frequency. This problem is analogous to 
that posed by a non-orthogonal Hamiltonian in the case 
of electrons, if we replace the energy, Hamiltonian, and 
overlap, by the square frequency, dynamical matrix, and 
mass matrix, respectivelji2ii2^. 

The presence of matrix M instead of the identity ma- 
trix introduces mathematical difhculties, if Green func- 
tions are used. The issue has been elegantly dealt with in 
the theory of non-orthogonal tight binding^^-: the "non- 
orthogonal" problem is equivalent to the "orthogonal" 



(c^2i + M-1/2KM-1/2)u = 



(6) 



The elements of the " orthogonalized" dynamical matrix 
elements are explicitly defined, as 



(7) 



and the dynamic equation is recast as 

{uj^I-K)u = 0, (8) 
which now can be treated using Green function methods. 

2. Local heat current density 

We are faced with the problem of calculating the local 
heat flow. The following paragraphs derive the necessary 
equations. The section following this one will express all 
the equations in terms of Green functions, so that there 
is no need to calculate individual wave functions. 

Let us consider an arbitrary wave (not necessarily an 
eigensolution) Ui(t), propagating in the system. The to- 
tal energy can be expressed as a sum over each degree of 
freedom, E — J^i ^i, with 
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Using Miili — —kijUj, the local change of energy with 
time is 
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The local current between each pair of local degrees of 
freedom is thus naturally defined, as 



(11) 



For a given phonon of frequency u we rewrite u and ii 
in terms of the complex wave 4>{t) = -06*"*: 



u,{t) = Re[Mt)]/^^ = (^f{t)/^/I'U (12) 
U^{t) = -LJIm[<j>,{t)]/VM^ = -oJ(f>i (t)/VM^ (13) 

Hence, the current associated to that particular phonon 
between the i and j local degrees of freedom is 



UJ 



-Im{(t)*hj(t)j) 



(14) 



The (j) are solutions of the eigenvalue problem, Eq. ((SJ. 
The normalization condition for the phonon amplitude 
follows from equating the wave's energy to fiuj^ as 

^|0,P = 2n/c^ (15) 



3. The total current 

We now proceed to derive the expression for the total 
heat current given exclusively in terms of Green func- 
tions. To this end, we subdivide the whole system into 
three regions, as indicated in Fig. ^ The two interfaces 
define four special groups of atoms labelled a, a, b and 
(3, as depicted in the figure. An atom belongs to one of 
these groups if there is a non-zero element of the dynam- 
ical matrix linking it to another atom on the opposite 
side of the interface. In what follows, a bar on top of a 
symbol, say Tp, indicates the one column matrix formed 
by the values of Lpi at each degree of freedom of the par- 
ticular group of atoms considered. A parenthesis denotes 
the subsystem: for example ^(a) is formed by values of 
Lf at the atoms of subgroup a only, and its dimension is 
3 times the number of atoms in a. Similarly, a bar below 
the symbol, say (/j, denotes a one row matrix. We denote 
the traveling wave solutions of the total system by 
The total current in terms of them is then given by (see 
Eq.llil) 



(16) 



where is the number of phonons in state n. 

Using standard scattering theorjsi^ one can express the 
eigenfunctions of the total system, <&„ in terms of the re- 
tarded Green function matrix of the total system, G, and 
the eigenfunctions of the decoupled system at the left of 
interface 1 (that is, with the coupling between a and a 
set to zero) . Denoting these wavefunctions for the decou- 
pled system as 0(a), one has for the waves propagating 
from a towards /3, 

(a) (17) 

^n(/3) = G/3akaQ0„(a) = g/S/sk/JfeGfaakaQ 0„ (a) , (18) 



4 





FIG. 1: Nomenclature used in Section III B 31 The semi- 
infinite leads 1 and 2 are projected onto subsets a and (3, 
and the total Green function for susbsystem 3 only needs to 
be computed between atoms of subsets a and b. 



where gaa and gfjp are the retarded Green function ma- 
trices corresponding to the decoupled systems (i.e. with 

kaa = and kb/i = 0). 

Inserting these expressions in the equation for the cur- 
rent we have 

J = 

(a) 

71 

= ATV^Wn X 

n 

xTr [^/TO{^„(a)^(a)kaaG*j,kb/3g/3^k^bGf,akaQ} , 

(19) 

where the sum extends to all states of the uncoupled 

subsystem 1, and 
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is the occupation difference between phonons travelling 
forwards and backwards. 

It is not feasible to explicitly compute the eigen- 
functions 0„ for a semi-infinite system. This can be 
avoided by expressing everything in terms of Green func- 
tions, which allow to project semi-infinite systems on the 
atomic groups considered. From the normalization con- 
dition, Eq. (|15l) . and the well known relation between the 
green function and the density of statesi^, it follows that, 
at any infinitesimal frequency square interval, {e, e-l-de}, 

w„0„(a)0* (a) = 2-/m[gQa(e)]de, (21) 

Z-^ TT 

where t = ijp. Using this relation, the expression for the 
current, Eq. H19() . becomes 
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where the transmission function S((jj) has been defined. 
The thermal conductance, Eq. ||2J), follows straightfor- 
wardly from Eq. H22|l . Now we only need to compute the 
Green functions g and G, and then use Eqs. (|23|) and 
(I22|l to calculate thermal currents. 



4- Calculation of the Green functions for infinite systems 
without translational periodicity 

The system we are treating is infinitely extended in the 
heat propagation direction, but it does not have transla- 
tional symmetry. Therefore, the calculation of the Green 
functions requires the use of projection techniques. The 
most efficient way to obtain the Green function of the 
semi-infinite systems 1 and 2 projected at subsystems 
a and /3 is to use the decimation techniqueS&, which is 
based in a renormalization procedure. After the gaa and 
gpp are calculated in this way, one can calculate the to- 
tal retarded Green function of the system everywhere in 
subsystem 3 (Fig. ^ , as 



G(cj2) = [uH - K - Si(cj2) - ^2(^2 



(24) 



where the self-energy matrix Si is defined as ^aa^aa^aa 
on the elements belonging to subset a, and zero other- 
wise, and S2 is defined similarly on subset b. 



C. Self-consistent flux model for efficient 
computation of boundary scattering 

If we know the dynamical matrix, the formalism ex- 
plained in the previous paragraphs is capable of account- 
ing for any type of scattering mechanism other than an- 
harmonicity (for the role played by anharmonicity, see 
appendix). It is straighforward to include isotopes by 
varying the mass of given atoms; substitutional impuri- 
ties may be introduced by varying both the mass of se- 
lected atoms and the links to their neighbors; vacancies 
would be included by removing single atoms and allow- 
ing the system to relax; etc. We might thus be tempted 
to treat boundary scattering in this same spirit, by in- 
cluding surface adsorbates or imperfections that can lead 
to effective boundary scattering in the nanowire. Adding 
single adsorbed impurities is a rather simple matter, but, 
could one consider an amorphous overlayer on the same 
footing? For example, it is known that Si nanowires are 
covered by an amorphous Si02 overlayer, several nm in 
thickness^, and it is important to know how this over- 
layer affects the phonon scattering. A brute force ap- 
proach would be to input the atomic positions of all con- 
stituents in the silica overlayer, as well as the wire, and 
then calculate the dynamical matrix of the whole system 
and compute its transmission. 

This procedure, however, is a nearly impossible task in 
practice. The added overlayer considerably increases the 
size of the system, and it also results in an extremely ir- 
regular transmission function, which has to be computed 
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FIG. 2: Left: scheme of a linear chain system showing the 
nomenclature used in the self consistent flux technique of Sub- 
section lll Ul Right: scheme of 1-D wire with Bethe lattices of 
coordination 4 attached. 



at a very large number of different frequencies. The in- 
tegration of such a transmission in Eq. (j^J , consequently, 
becomes cumbersome. To solve these difficulties, we have 
developed the self-consistent flux model. This model's 
basic idea consists in taking the limit of an infinitely 
thick overlayer. In the same spirit of the Bethe lattice ap- 
proach, as used for amorphous materialsi^, we consider 
the nanowire connected to Bethe lattice branches of the 
amorphous material (see Fig.[2J). The phonon mean free 
path in amorphous materials is of the order of the in- 
teratomic spacingiSi. Therefore, the lack of connectivity 
in the Bethe lattice approach does not represent a prob- 
lem. The adequacy of the Bethe lattice approach for 
the study of amorphous systems has been investigated in 
refs. l2(ill2l 

When the overlayer is substituted by an infinitely thick 
one, an essential difference appears. Now, the self- 
energies associated to the overlayer do not correspond 
to clusters. Therefore, they no longer are purely real, 
but contain a finite continuous imaginary part. In other 
words, if one calculates the transmission of the system by 
Eq. (|23|l . part of the phonon flux will be effectively lost 
through the side branches, since now they are inflnite and 
allow for the current to go away without returning. How- 
ever, this is not physically correct, for we know that the 
overlayer is finite and the heat current must eventually 
come back to the wire and escape through the contacts. 
In a steady state, this means that we must allow for an 
equivalent ammount of phonon flux to re-enter the sys- 
tem through each of the side branches. One may flnd 
a self-consistent solution to this problem in the form of 



an iterative procedure. This is indeed not necessary, for 
the self-consistent solution can be directly found by an 
algebraic method, as we show below. 

To clarify the nomenclature in the following discussion, 
a general scheme of the nanowire with part of its length 
coated is shown in Fig. 12 left. The represented nanowire 
is three dimensional, thus each black dot in the figure 
stands for a group of atoms, and each line joining dots 
stands for a set of interatomic links. The dotted line en- 
closes the wire's groups of atoms that are coated. Bethe 
lattices are attached to those groups of atoms, as shown 
in the figure. Phonons enter the system through group 1, 
and exit at the other lead through group n. First of all, 
we obtain the partial transmissivities, T, between every 
pair of branches in the system, using Eq. (|23|l . 

4 - _ - 

Tij = —Tr[ku'Im[g^,^,]k,,^G*jkjj,Im[gj,j,]kj,jGji]{25) 

where i and i' label neighboring atomic subsets belonging 
to the branch and wire respectively, for the first branch of 
the pair, and equivalently j and j' for the second branch. 
N is the number of degrees of freedom in node i' or j' . 
(All nodes are assumed to contain the same number of 
atoms for simplicity.) 

We recall that each link to a branch in Fig. [3 may 
comprise interactions with several atoms, in the three 
geometrical directions, so the Green functions above are 
matrices. The transmissivity from one branch to itself 
is not defined by the above formula. Instead, we must 
define it as the reflectivity, 

Tu^l-Y.'^^i- (26) 

We deflne the "injected" and "outflowing" flux arrays, 
ji and jo as the fluxes coming "into" and "out of the 
system via its surrounding nodes (see Fig. |2J). The total 
flux array, j, is then 

J=J^-Jo (27) 

and it must fulflU 

j(l)=Ji(l)-Jo(l)=Jo(") (28) 

because all the flux coming in through node 1 must flnally 
exit through node n. Also, there is no flux injection from 
node n, so 

Un) = (29) 
All the other modes must fulfill 

j{l)^Ml)-Jo{l) = 0,Vl^l,n (30) 
We can express Eqs. (|28I29I30|I in matrix form as 

= Ajo (31) 
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with 



A = 
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(32) 



On the other hand, the outflowing currents are related 
to the injected ones by 



which implies 



jo = ^Ji 



jo = TAjo 



(33) 



(34) 



Therefore, the self-consistent array jo is proportional 
to the eigenvector vi of TA with eigenvalue 1, such that 
TAvi = vi. The existence of eigenvalue 1 is guaranteed 
by the way we have defined Tu (Eq. H26|l) and by the fact 
that columns 1 and n in matrix TA coincide. 

The transmission corresponds to the outgoing flux at 
node n when there is one unit of flux inciding per de- 
gree of freedom in node 1. Thus, imposing ji{l) = N = 
jo(l) -|- join), where the second equality follows from 
Eq. |(2H1), gives 



N 



Jo 



vi{l) +vi{n) 
and the self-consistent transmission is 



(35) 



S = Join) = N 



vi{n) 



vi{l) + vi{n) 



(36) 



It can be easily verified that for the uncoated case the 
calculation using Eq. H36I) coincides with that of Eq. (|23|l . 

In the next section we explicitly demonstrate that the 
SCFlux approach yields the same result that the calcula- 
tion using large clusters attached to the wire. Therefore, 
the SCFlux technique is of tremendous importance, since 
it enables to study the otherwise untractable problem of 
thermal transport through thick-amorphous-coated real 
wires. 



III. RESULTS 

New physics regarding how the amorphous coating in- 
fluences the thermal transport through a wire can now be 
learnt using the above method: Does the coating affect all 
frequencies equally, or are there important differences in 
the scattering as a function of frequency? Does diffusive 
transport arise? And, if so, how long a segment has to 
be coated in order for transport to diffusive rather than 
ballistic? How does transport depend on the coupling at 
the wire-coating interface? Etc. 



In the next two subsections we present results for two 
different systems. First, a one dimensional model, which 
is simple enough to allow for a thorough investigation 
without encountering size limitations, and at the same 
time displays many general features that apply to real 
systems. Afterwards, we present results for real 3-D sili- 
con nanowires. 



A. Study of a one dimensional model 

A wide range of interesting phenomena is obtained al- 
ready in the case where there is only one degree of free- 
dom per atom, in a linear chain. Bethe lattices with also 
one degree of freedom per atom are attached to a section 
of the chain. The wire's properties are determined by the 
value of a spring constant Kw- For the Bethe lattices, 
we assume them to be linear chains (i.e. Z = 2 in the 
nomenclature of ref . il&l , and they are characterized by 
their spring constant Kg. The third constant involved is 
the coating- wire interaction, Kbw, connecting the end 
of the Bethe lattice to the wire. 



1. Comparison between SCFlux and cluster calculations 

First of all, we provide a particular example that 
demonstrates the equivalence between the SCFlux calcu- 
lation, and the much more computationally demanding 
calculation using large finite clusters for the coating. The 
transmission was computed for a system with ten clus- 
ters attached to it at ten consecutive atoms. Each of the 
clusters is composed of sixty overlayer atoms. The trans- 
mission, calculated via Eq. H23|) . is shown in Fig. |2Ia). 
We shall abbreviately refer to this curve as the "clus- 
ter transmission" . A large number of closely spaced an- 
tiresonances occurs, due to the finite size of the clusters. 
Thus, the curve had to be computed at many different 
frequency points 42i 

The transmission given by the SCFlux method is 
shown as the thicker line in Fig. Olb) . We shall refer to 
it as the " SCFlux transmission" . Its smooth shape has 
little visual resemblance with the cluster transmission. 

We can now compute the thermal conductance as 
a function of temperature, using both the cluster and 
SCFlux transmissions, and compare the results. We de- 
fine the dimensionless thermal conductance as 
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where lju = 2Kb is the Debye frequency of the infinite 
1-D atomic chain. Using Eq. lO, the dimensionless con- 
ductance can be evaluated in terms of the dimensionless 
temperature T = j^j^T, and the dimensionless frequency 
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(38) 
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FIG. 3: Comparison between the transmission calculated by 
using finite clusters as scatterers or the SCFlux method with 
infinite Bethe lattices, a) Cluster transmission. (10 clus- 
ters comprising 60 atoms each are attached to the chain.) 
b) Thick line: SCFlux transmission; thin line: the curve in 
"a)" averaged in a range of frequencies around each point. 
As more atoms are added to the clusters, the averaged curve 
approaches the SCFlux result more closely. 



FIG. 4: Comparison between the dimensionless thermal con- 
ductance as a function of the dimensionless temperature, cal- 
culated from the cluster transmission (dashed line), and from 
the SCF/iis transmission (solid line). Although the two trans- 
missions have very different shapes (c.f.Fig. I^J, the conduc- 
tances obtained from them are nearly the same. The conduc- 
tance obtained from a constant transmission function is also 
shown (dotted line) to stress the importat role played by the 
frequency dependence of the transmission function. 



The dimensionless conductance calculated from the 
cluster transmission is shown as the thick dashed line 
in Fig. 21 as a function of the dimensionless temperature. 
Now this curve can be compared with the thermal con- 
ductance calculated using the SCFlux transmission. Re- 
markably, the latter yields the thick solid curve in Fig. 0) 
almost exactly the same thermal conductance curve than 
that computed from the cluster transmission. The rea- 
son why the two apparently very different transmissions 
yield virtually the same conductance can be better un- 
derstood by averaging the cluster transmission in a range 
of frequencies, so as to smooth out its structure. We have 
convolved the curve in Fig. OJa) with a truncated Gaus- 
sian kernel extended ±wd/60 around each point. The 
cluster transmission averaged in this way is shown as the 
thin line in Fig. |31b), and is very close to the SCFlux 
transmission. 

Thus we see that the transmission calculated with the 
SCFlux method does indeed yield the same thermal con- 
ductance as if we include the thick coating in the form 
of large finite clusters. But the SCFlux method requires 
a much smaller computation, because it yields a smooth 
transmission curve, while the cluster calculation needs to 
evaluate a far larger number of points in order to account 
for the closely spaced antiresonant structure associated 
to the finite clusters. 

To quantify the difference between the real conduc- 
tance and the one obtained using a frequency independent 
boundary scattering, the latter is also plotted in Fig. ^ 
as the thin dotted curve. The value of the constant trans- 
mission has been fixed to yield the same limiting conduc- 
tance in the lu oo limit as the real conductance curves 



(see inset of Fig. 0J. As the figure shows, important dif- 
ferences between the oversimplified constant transmission 
model and the real transmission calculation exist up to 
well above the debye temperature. 

2. Study of the transmission function 

Now we proceed to study the behavior of the trans- 
mission as a function of the different physical parame- 
ters in play. More concretely, we want to know how the 
transmission depends on frequency lo, coated length L, 
and strength of the interface spring constant Kbw- The 
analysis will show us how the system goes from ballis- 
tic to diffusive behavior, as a function of these physical 
parameters. 

We set Kw to be the unit of spring strength. The 
mass of all atom s, M, is defined as 1 also. Frequencies 
are in units of Kw /M. The coated length L is given 
in lattice units, being equivalent to the number of Bethe 
lattices attached. 

First let us take a look at the behavior of the transmis- 
sion function with frequency, in Fig. El (for Kw/Kb = 1 
and Kbw = !)• We notice that the effect of boundary 
scattering is to reduce the transmission rather uniformly 
throughout the whole frequency range. However, it is 
apparent that the usual assumption of a constant scat- 
tering ratei2iSL22iSS is only a rough approximation to the 
real scattering in nanowires. The low frequency limit 
behavior displays an important feature: for wires where 
the coated segment is long enough (n > 3 in the case 
shown), the transmission near w ~ decreases with in- 
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FIG. 5: SCFlux transmission as a function of frequency for 
infinite 1-D chain, for coated segments of different lengths, L 
(= number of lattices attached). 



creasing frequency. This appears to be the general case, 
as we will see in the next section. It has a direct physical 
consequence: the slope of thermal conductance at very 
low temperatures decreases below its T = value when 
temperature increases. This fact has been shown to hap- 
pen also in the case of surface roughness at non-coated 
wiresii (the well known "dip" of refj2r) We explicitly 
show this for Si nanowircs in the next subsection. 

In Fig. [SI we have seen that the shape of the transmis- 
sion function varies as a function of the coated length. 
Does it attain a limiting shape? And, how is this limit 
reached? This is the subject of the diffusive-ballistic tran- 
sition. The way it takes place will become clear from 
the study of the transmission as a function of coated 
length, L. In principle, it is not obvious what type of 
length dependence one could expect. For pure disorder, 
for example, an exponential decrease in transmission is 
obtained in 1-D systems due to localization'^^. As an- 
other example, for wires with rough edges, an anomalous 
L~^/^ thermal conductance has also been reported un- 
der certain conditions*. For amorphous-coated wires, the 
asymptotic behavior turns out to be ^ L^^, as we now 
show. 

Analogously to ref. l3lL we write the transmission as 

S = So(l + L/A(L))-i (39) 

This is merely a definition of A(L), and it does not impose 
any limitations on the form of the transmission. For the 
present 1-D system, the uncoated wire's transmission, 
5o, is equal to 1 at all frequencies. Now, the form of 
A(L), plotted in Fig. |H| tells us the length dependence 
of the transmission. New physical insight can be gained 
from this figure. First, we see that X{L) quickly tends to 
an asymptotic value. This means that, for amorphously 
coated wires, the long L behavior of the transmission 
follows an L^^ dependence^. More interestingly, the 
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FIG. 6: Plot of A(L) = (Hq/H - 1)-^L as a function of the 
coated length L, for different frequencies. X{L) asymptotically 
converges to a limiting relaxation length value, Aoo((jj) for long 
enough L. 

length it takes X{L) to saturate to its asymptotic value 
becomes longer as the frequency gets closer to zero. Phys- 
ically, this implies that the lowest frequency phonons 
can still travel nearly ballistically for system lengths at 
which higher frequency phonons are transported in an 
almost totally diffusive fashion. In a practical measure- 
ment of nanowire thermal conductivity, this might lead 
to a wrong estimation of the low temperature thermal 
conductivity, if the wire were not long enough. 

In the diffusive regime, we can define a frequency de- 
pendent relaxation length, Aoo, as the limit 

Xoc{io) = limL^ooK^)- (40) 

A very interesting phenomenon arises when we study 
how this relaxation length depends on the strength of 
the spring constants coupling the wire to the amorphous 
coating. We numerically computed Aoo at frequency 
Lu = 0.05a;/), as a function of the wire-coating coupling 
strength. The result is plotted in Fig[7| For decreasing 
values of the wire-coating coupling, Kbwi the relaxation 
length quickly increases, as expected. However, if the 
coupling is made stronger, Xoo reaches a minimum and 
then it increases again as the coupling is increased, ap- 
proaching an asymptotic value. The existence of this 
minimum relaxation length can be interpretted in micro- 
scopic terms: there is an optimum value of the coupling 
that maximizes the heat-flow exchange to and from the 
coating, thus maximizing the scattering. For very weak 
couplings the heat-flow exchange reduction is intuitively 
obvious. For very strong couplings on the other hand, 
the linked atoms would behave more like a rigid cluster, 
thus acting as a hard wall that also confines the phonons 
to the wire rather than let them enter the amorphous 
region. From a macroscopic point of view, the relaxation 
length is minimum when the specularity factor is 0, in 
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FIG. 7: Phonon relaxation length Aoo as a function of the 
coupling between the wire and its coating lattices. It reaches 
a minimum and then increases, asymptotically approaching 
a limiting value. (In the particular case shown here, lu — 
O.OSojD, and A 2.7.) 
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FIG. 8: Phonon dispersion relations calculated for bulk Si 
using Harrison's potential. 

the Casimir limit"^^. The minimum of Aoo obtained thus 
translates to this hmit, and specularity increases in both 
ways around this point as a function of Kbw- 



B. Results for Si nanowires 

The advantage of the atomistic description used in this 
paper is that we can study wires made of " real" materi- 
als. In other words, we are dealing with the real phonon 
spectrum of the system: the bands are dispersive, and 
their frequency extension is finite, in contrast with the 
infinitely extended spectrum in a continuum elastic sys- 
tem. 

We proceed to study the transmission of phonons in 
real Si nanowires, with an amorphous material coating 
the wire. The phonon dispersion relations of Si are fairly 
well reproduced with an interatomic potential that in- 
cludes only two and three body terms. We use Harrison's 
potentials. The adequacy of this potential for Si is as- 



FIG. 9; Cross sections of the three Si nanowires studied: 
"2x2", "3x3" and "4x4". 

sesed by the satisfactory dispersion relations calculated 
with it. We show them for bulk Si in Fig. |S1 Despite 
having only two parameters, the shape and position of 
the bands are in reasonably good agreement with exper- 
imental results'^'*-^^. 

The two body contribution to the energy is 

5Eoi^,J)^lCo^-^^^^^^ (41) 

for every pair of nearest neighbors i and j, where dij 
is the distance between the atoms and do is the lattice 
equilibrium distance. The three body contribution is 

5E,izJ,k) = ^C,Sel^, (42) 

for every pair of bonds joining atoms i, j and fc, where 
SQ is the deviation with respect to the equilibrium angle 
between the two bonds in the lattice. Constants Cq = 
49. ley and Ci = l.Q7eV are taken from table 9-1 of 
ref.isa 

For the coating we again consider Bethe lattices con- 
nected to each vibrational coordinate of the boundary 
atoms. The parameters of the Bethe lattice are chosen 
to reproduce the frequency range of the known d.o.s. of 
silica™, i.e. 2 x 10^^ Hz. A coordination number Z = 4 
is assumed. We have considered (Oil) wires of three dif- 
ferent widths, shown in Fig. |3 The wires are coated 
by attaching Bethe lattices only to the outermost atoms. 
For the 2x2 wire all atoms are at the surface. The 3x3 
wire has 16 surface atoms and 2 core atoms. The 4x4 
wire has 24 surface atoms and 8 core atoms. The phonon 
dispersion relations for each of them is shown in Fig. 1101 

The procedure for the calculation is as follows. First, 
the nanowire lattice is defined, and the dynamical matrix 
is obtained for the periodic unit cell, and for the link be- 
tween two neighboring unit cells. The Bethe lattices are 
included as a self-energy projected onto their last atom. 
The self-energy due to the uncoated parts of the wire is 
calculated in the way of ref. |2^ Then the Green func- 
tion of the system is calculated (Eq. (|24|l ). The partial 
transmissivity between each pair of branches (Eq. (|25|l ) is 
then calculated, and the SCFlux transmission (Eq. (|36|) ) 
is computed. 



10 




D 0.2 0.4 0.6 D.8 1 

k [n/a] 




0.2 0.4 0.6 0.8 1 

k [Kh] 



0.8 



3 








0.2 D.4 0.6 0.8 1 



k [Tz/a] 

FIG. 10: Phonon dispersion relations for the three Si 
nanowires in Fig. |5| 



We show how the phonon transmission depends on fre- 
quency, for different values of the coated length. Fig. ITTl 
shows as the coated length increases from to 5 

coated unit cells. In these curves, all three spatial direc- 
tions of each of the surface atoms are attached to Bethe 
lattices. The qualitative behavior is an overall uniform 
decrease of the transmission curve, that maintains the 
peaks and curve features of the uncoated wire. This is 
in agreement with the usual assumption of a constant 
scattering rate. However, quantitatively the boundary 
scattering rate is not frequency independent, as it is ap- 
parent from the relative heights of the peaks at different 
coating lengths. In the absence of coating, the transmis- 
sion curve at the lower end of the frequency spectrum has 
a value of 4. This is always the case, since there are four 
lowest frequency branches at all nanowires: one dilata- 
tional, one torsional and two flexural branches^!. The 
two flexural branches have a quadratic rather than linear 
dispersion, as we can see in Fig. ()10|l . being a consequence 
of the soft character of those modes. 

The transmission at the lowest frequency decreases 
more slowly than that around the first peak {^-^ 20THz), 
as the coated length becomes larger. Thus scattering 



is weaker for the lowest lying modes than for the bulk 
of transverse modes. This results in a flattening of the 
transmission as L increases. The narrower the wire, the 
more pronounced this effect becomes. 

The low frequency behavior of the transmission func- 
tion also has important physical implications. Similarly 
as what we saw in the 1-D example of Section IIII Al . 
for real Si nanowires the transmission function also de- 
creases below its = value for small values of the 
frequency, and increases again as the frequency rises fur- 
ther (see inset in Fig. I12|l . The result is that, at low, 
flnite temperatures, the thermal conductance divided by 
T decreases with respect to its zero temperature value, 
thus displaying a dip (see Fig. ^J. A dip in a{T)/T 
has been reported in ref. 0, and theoretically explained 
in ref. ^3 to be the result of surface roughness scattering. 
In our case of amorphous-coated nanowires we see that 
an analogous dip occurs, although the system considered 
here is basically different from that one, and the con- 
ductance reduction phenomenon cannot be considered to 
be the same one. Some essential differences exist in the 
case of scattering by amorphous coating with repect to 
the roughness scattering case. One is that, at zero fre- 
quency, the scattering rate is zero in the case of roughness 
scattering, while it is different than zero in the case of 
amorphous-coating scattering. As a result, it is possible 
to measure a conductance close to 4 times the quantum 
of conductance at uncoated wires, but for wires with a 
thick amorphous coating the measurement will in general 
yield smaller values at the very low temperature limit. 
Another difference pointed out earlier is the exponential 
decrease of conductance with length in the rough wires, 
as compared with the inversely proportional dependence 
in the long length limit for amorphous-coated wires. 

We have also calculated the relaxation length, Aoo 
(Eq. if^ l. as a function of the Bethe lattice coverage, 
in Fig. El The abscissa corresponds to the number of 
Bethe lattices attached to each pair of surface atoms of 
the 4x4 nanowire. Since there are three spatial direc- 
tions, we can attach up to 6 lattices to each of the pairs. 
We show Aoo as a function of the number of directions at- 
tached, for different ways of attachment. The directions 
in which the Bethe lattices are attached, and the number 
of them, are depicted schematically near the points in the 
graph. The calculation shows that the relaxation length 
is strongly dependent on the character of the bonding to 
the coating material. For a densely covered wire, the re- 
laxation length becomes of the order of the wire's core 
diameter. On the other hand, if the surface bonds to the 
coating are less dense, with one dominant direction and 
two much weaker links in the other two directions, then 
the relaxation length can attain values between 2 and 4 
times the diameter. When there is only partial bonding of 
part of the surface atoms, the relaxation length increases 
further. The way in which A approaches its saturation 
value, Aoo is shown in the inset of Fig. E| With only 
one link per surface atom (solid curve), it is necessary to 
coat a considerable length in order to attain the diffu- 
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FIG. 11: SCFlux transmission for Si nanowires. Left: transmission of the 2x2 wire for 0, 1, 2, 3, 4, and 5 unit cells coated. 
Center: transmission for the 3x3 wire with 0, 1, and 5 unit cells coated. Right: same, for the 4x4 wire. 
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FIG. 12: Low temperature behavior of the thermal conduc- 
tance divided by temperature, for the 4x4 nanowire with 5 
unit cells coated, showing a dip similar to the one studied 
in im (The ordinate axis unit is ^Ti/k%.) Inset: detail of 
the transmission function 'S.{lu) at low frequencies, for the 4x4 
wire with 1 (dotted line) and 5 (solid line) unit cells coated. 



sive limit. With more interface links, the diffusive limit 
is reached at shorter coated lengths. It is observed that, 
not only the number of links, but also their orientation, 
has an influence in the relaxation length. 

A larger number of interface links does not necessarily 
imply a shorter relaxation length. A fully covered surface 
with three links per atom, for example, displays a larger 



relaxation length than a surface with only half of the 
surface atoms directly attached to the coating (Fig. I13|l . 
This is another manifestation of the phenomenon noted 
in Subsection IIII Al where we saw that a very strong in- 
teraction between wire and coating can result in longer 
relaxation lengths. 



IV. CONCLUSIONS 

We have presented a formalism to calculate the phonon 
transmission function of an atomically described system. 
With this formalism we have studied the problem of 
phonon transport in nanowires with part of their length 
coated by an amorphous material. A self-consistent ap- 
proach has been developed that uses Bethe lattices to 
treat the amorphous material. The efficiency of the ap- 
proach, and its equivalence to a much more computation- 
ally demanding cluster calculation, have been explicitly 
illustrated via direct comparison, by an example. 

We have studied two examples of amorphous-coated 
wires: a 1-D model, and silica coated Si nanowires. In 
the examples studied, specific physical phenomena for 
coated systems have arosen: 

1) Ballistic- diffusive transition: We explicitly showed 
the evolution of the phonon transmission in the transi- 
tional regime, before A(a;,L) asymptotically approaches 
its limiting value. In the limit of long coated length the 
transmission attains a 1/L behavior. (In contrast, for 
rough-boundary uncoated wires, an exponential decrease 
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FIG. 13: Phonon relaxation length Aoo in 4x4 wire at fre- 
quency uj — 16 THz, for different ways of attachment between 
the coating and the wire surface atoms, in units of the wire's 
diameter D. The scheme near each datapoint represents the 
directions in which the Bethe branches are attached to each 
surface pair of atoms. The inset shows A versus the coated 
length L (in unit cell lengths), for three particular attachment 
configurations. 



with length takes placaii.) 

2) Low temperature conductance: near w = 0, the 
transmission is found to decrease with respect to its 
UJ — value, as the frequency is increased. As a result, 
amorphous-coated wires display a low temperature dip in 
the plot of the thermal conductance divided by tempera- 
ture, having certain similarities with the dip reported for 
rough uncoated wires^^. 

3) Coating-wire interface structure effects: a new phe- 
nomenon is found, which is specific of coated wires: as 
the interatomic bonds between wire and coating atoms 
are made stronger, starting from zero, the transmission 
of a long coated wire (or equivalently, A(L, u)) decreases, 
but only up to a certain limit; once the bond strength 
reaches a certain threshold, the transmission increases 
again as the bond is made stronger. When the bond 
strength tends to infinity the relaxation length asymptot- 
ically saturates to a finite value. An analogous behavior 
was found to take place as a function of the number of 
interface links. The structure of the interface was shown 
to have an important effect on the relaxation length. 

It is hoped that the method and results presented 
here will help to estimulate further atomistic research 



in the field of phonon transmission through nanowires, 
as well as experimental research on amorphous-coated 
nanowires. 
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APPENDIX: THE ROLE OF ANHARMONICITY 

In one dimension, there is no three phonon anharmonic 
scattering, since it is not possible to have a process in 
which all three phonons belong to the same polarization 
branch^. Therefore, anharmonicity is not an issue in the 
first example, given in section Fill Al 

The role of anharmonicity in the thermal conductiv- 
ity of Si nanowires has been studied in ref. For the 
wires studied in section UlI Bl its influence is negligible. 
An estimation can be made using the Mathiessen rule. 
The total relaxation length, expressed in terms of the 
boundary and anharmonic relaxation lengths, is: 



A = l/(A-i + A-i). 



(A.l) 



Now, Afe ~ Z?, where D is the system's thickness (around 

2 nm, for the wires considered here). Also, Aa ~ 
c/(1.73 X 10-i^;|rw2e-i37.3K/T-)^ .^i^gj.g ^ jg ^j^g gi gpgg^ 

of sound'^^. The role of anharmonicity is largest for the 
highest frequencies. But even for the highest frequency 
available in Si, ^ 95 THz, the anharmonic relaxation 
length at room temperature is still Aa ~ 20 nm. This, in 
the case of the nanowires of section lTlI Bl using eqn. I|A.1|) 
and the estimate just given for Ab, affects the transmis- 
sion by less than 10%. (Lower frequencies are much less 
affected.) The influence on the thermal conductivity is 
much smaller than this, since high modes contribute sig- 
nificantly less than the lower energy modes. Experimen- 
tally it has been observed that even for nanowires 37 nm 
thick, the effect of anharmonicity on the thermal conduc- 
tivity is minimal compared to that of boundary scatter- 
ing, up to 300K. This is evident from the fact that the 
thermal conductivity stays constant with T near room 
temperature, rather than decrease (see refs. I3&.39) . 
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